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In their constructions of system of quantum stochastic differential equations, mathemati- 
cians and / or several physicists interpret that the function of random force operator is to 
preserve the canonical commutation relation in time, i.e., to secure the unitarity of time 
evolution generator even for dissipative systems. If this is the case, it means physically 
that the origin of dissipation is attributed to quantum non-commutativity (quantum- 
ness). The mechanism that the mathematician's approaches rest on will be investigated 
from the unified view point of Non-Equilibrium Thermo Field Dynamics (NETFD) which 
is a canonical operator formalism of quantum systems in far-from-equilibrium state in- 
cluding the system of quantum stochastic equations. 



1 Introduction 

There are arguments EHHHiD that the function of random force operator is to 
preserve the canonical commutation relation in time. The contents of issue are 
the following. The time-evolution of free Boson operators is given by da{t)/dt = 
(t), da*{t)/dt = iu;a*(t). The canonical commutation relation [a, a*] = 1, at 
time t = preserves in time, i.e., [a(t), a\t)] = 1. If a relaxation is introduced 
simply by 

da(t)/dt = -iua(t)-Ka(t), da f (t) / dt = iuja f (t) - na f (t) , (1) 

the canonical commutation relation decays as [a(t), a*(t)] — e~ 2fit . This incon- 
venience is secured by introducing random operators F(t) and F*(t) which are 
assumed to satisfy [a, F^(t)] = 0, [aJ, F(t)] = 0, etc. for t > 0. The solutions of 
Langevin equations 

da(t) /dt = -iua(t) - na(t) + F(t) , da) (t) /dt = iuja) (t) - no) (t) + F f (t) , (2) 

are given by a(t) = ae~ tuJt - Kt + /* dt'F(t')e- iu '^- t ^ +K( - t+t '\ al(t) = a ^e luJt - Kt + 
L dt 1 F^ '(i^e*^* -4 )+ K (*+* ). Then, one knows that the canonical commutation rela- 
tion preserves in time if the condition 

l = [a(t), at(i)]=e- 2K *|l + ^ db x J dt 2 [F(t 1 ), p^t^e^-^+^A (3) 

is satisfied. It is realized by the commutation relation among the random force 
operators: 

[Fit), Ft(t')] = 2nS(t - t'). (4) 

■Wliese stochastic differential equations should be interpreted as those of the Stratonovich 
type Ii3, since we perform calculations as if they were ordinary differential equations. 
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The above argument seems to show us that the function of random force opera- 
tors is to recover the unitarity of time-evolution generator rather than to represent 
dissipative thermal effects. If it is correct, doesn't it mean, physically, that the ori- 
gin of dissipation can be attributed to quantum non-commutativity (quantumness)! 
There are, however, physical systems described by classical mechanics, or quantum 
systems with commutative random force operators. 

In this paper, withJjhe help Q£±he framework of Non-Equilibrium Thermo Field 
Dynamics (NETFD) □ffl0OLlL3'0, we will investigate the above argument in a 
systematic manner by means of martingale operator paying attention to the non- 
commutativity among random force operators. In section ||, the structure of the 
system of stochastic differential equations in classical mechanics is reviewed. In 
section |], the system of quantum stochastic differential equations within NETFD 
is introduced. In section || Boson system will be treated, which has a linear dissi- 
pative coupling with environment system within the rotating wave approximation. 
The mathematician's arguments will be studied from the unified viewpoint based 
on the canonical operator formalism of NETFD by changing the intensity of non- 
commutativity parameter A among a martingale operator. In section Boson 
system having a linear dissipative interaction between environment without the ro- 
tating wave approximation will be investigated with the help of NETFD. This is 
the case where the system has commutative random force operators. Section ^| will 
be devoted to some remarks. 

2 System of Stochastic Differential Equations in Classical Mechanics 

2.1 Stochastic Liouville Equation 

We will show the structure of the system of classical stochastic differential equa- 
tions E-3 starting with the stochastic Liouville equation 

df(u,t)=Sl(u,t)dtaf(u,t), (5) 

of the Stratonovich type with Q(u,t)dt — —{d/du)du where the flow du in the 
velocity space is defined by du — — judt + m~ l dR(t) . Here, the circle o represents 
the Stratonovich stochastic product E£l which is^defined in appendix A together 
with the definition of the Ito stochastic product E3. The increment of random force 
dR(t) is a Gaussian white stochastic process defined by the fluctuation-dissipation 
theorem of the second kind: 

(dR(t)} = 0, (dR{t)dR(t)) = 2mr/Tdt, (6) 

where 7 (> 0) is a relaxation constant, and T a temperature of the environment 
represented by the random force dR(t). Note that within the Stratonovich calculus 
dR(t) and f(u,t) are not stochastically independent: (dR(t) o f(u,t)) ^ 0. The 
average (• • •) is taken over all the possibility of the stochastic process {dR(t)}. 

By making use of the relation between the Ito and Stratonovich products (p2[), 
the stochastic Liouville equation (||) can be rewritten as the one of the Ito type: 

df(u,t) = (2(u,t)dtf(u,t), (7) 
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with n(u,t)dt = —(d/du)du where du is the flow in the Ito calculus given by 
du = —7 (u + m^ 1 T(d/du)) dt + m _1 <iR(i). Note that there appears temperature 
T in the flow, and that within the Ito calculus dR(t) and f(u,t) are stochastically 
independent: (dR(t)f(u,t))=0. 

The initial condition for the stochastic distribution function f(u,i) is given by 
/(it, 0) = P( u, 0) , where P(u,t) is the velocity distribution function defined below 
in subsection 2.3 . 

Note that the stochastic distribution function conserves its probability within 
the relevant velocity phase-space: J du f(u,t) = 1 



2.2 Langevin Equation 

The system described by the stochastic Liouville equation (|j) can be treated by the 
Langevin equation: 

du(t) = -ju(t)dt + m~ 1 dR(t). (8) 

This Stratonovich type stochastic differential equation does not contain the diffusion 
term, which is very much related to the way how physicists originally introduced 
the Langevin equation. 6 

It is worthwhile to note here that one could have introduced the Langevin 
equation within the Ito calculus of the form 

du(t) = -j[u(t) + mr l T{8/8u{t)))dt + rrT l dR{t), (9) 

which has a term with the functional derivative operator S/Su(t). In the system 
of quantum stochastic differential equations within NETFD, the Langevin equation 
of the Ito type, such as <M), can be introduced on the equal footing as the one of 
the Stratonovich type (H), although this was not the original motivation for the 
invention of NETFD (see section ||) . 



2.3 Fokker- Planck Equation 

In precise, the stochastic distribution function is given by f(u, t) = f(H(u, t), P(u, 0)) 
Taking the random average (•••), we have an ordinary velocity distribution function 
P(u,t) — (f(f2(u,i),P(u,0))) which satisfies the Fokker-Planck equation 

dP(u,t)/dt= (9/au)7[w + m~ 1 T(a/5u)]P(u,t). (10) 

This can be derived most conveniently from the stochastic Liouville equation (|7]) of 
the Ito type because of the orthogonal property mentioned in subsection [2.1[ 

The fluctuation-dissipation theorem (^) of the second kind is introduced in 
order that the stochastic Liouville equation f|]) and the Langevin equation (||) are 
consistent with the Fokker-Planck equation fllC|). 

'The Langevin equation was introduced by adding a random force term, such as m~ 1 dR(t) / dt, 
to a macroscopic phenomenological equation, for example, like du(t)/dt = —yu(t). 
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3 System of Quantum Stochastic Differential Equations 

3.1 N on- Equilibrium Thermo Field Dynamics 

In order to treat dissioatiy/? .mpituao. systems dynamically we constructed the 
framework of NETFD EMligOMJil It is a canonical operator formalism of 
quantum systems in far-from-equilibrium state which enables us to treat dissipative 
quantum systems by a method similar to the usual quantum field theory that ac- 
commodates the concept of the dual structure in the interpretation of nature, i.e. 
in terms of the operator algebra and the representation space. In NETFD, the time 
evolution of the vacuum is realized by a condensation of 7^7^-pairs into vacuum, 
and that the amount how many pairs are condensed is described by the one-particle 
distribution function n(t) whose time-dependence is given by corresponding kinetic 
equation (see appendix B). 

We further succeeded to construct a unified framework of the canonical operator 
formalism for quantum stochastic differential equations with the help of NETFD. 
To the author's knowledge, it was not realized, until the formalism of NETFD 
had been constructed, that one can put all the stochastic differential equations 
for quantum systems into a unified method of canonical operator formalism; the 
stochastic Liouville equation 113 and the Langevin equation within NETFD are, 
respectively, equivalent to the Schrodinger equation and the Heisenberg equation in 
quantum mechanics. These stochastic equations are consistent with the quantum 
master equation which can be derived by taking random average of the stochastic 
Liouville equation. 

3.2 Quantum Stochastic Liouville Equation 

Let us start the consideration with the stochastic Liouville equation of the Ito type: 

d\0 f (t)) = -iH u dt\0 f (t)). (11) 

The generator V f (t), defined by |0/(t)) = Vf{t)\0), satisfies dV f (t) = -iH f , t dt V f (t) 
with Vf(0) = 1. The stochastic hat-Hamiltonian Tif^dt is a tildian operator sat- 
isfying {iHf^dt)~ — iHf.tdt. Any operator A of NETFD is accompanied by its 
partner (tilde) operator A, which enables us treat non-equilibrium and dissipa- 
tive systems by the method similar to usual quantum mechanics and/or quan- 
tum field theory. Here, the tilde conjugation ~ is defined by (AiA 2 )~ = A1A2, 
(dAi + c 2 A 2 y = cjli + c* 2 A 2 , (i)~ = A, and (A^)~ = A* with A's and c's being 
operators and c-numbers, respectively. The thermal ket-vacuum is tilde invariant: 

|0/(<))~ = \o f (t)). 

From the knowledge of the stochastic integral, we know that the required form 
of the hat-Hamiltonian should be 

H f , t dt = Hdt + : dM t : (12) 

where H is given by H = Hs + in with Hs = H$ — Hs, and LT = flu + II d 
where LTr and Ub are, respectively, the relaxational and the diffusive parts of the 
damping operator 77. The martingale dMt is the term containing the operators 
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representing the quantum Brownian motion dB t , dB\ and their tilde conjugates, 
and satisfies (\dM t \) — 0. The symbol : dM t : indicates to take the normal ordering 
with respect to the annihilation and the creation operators both in the relevant and 



the irrelevant systems (see fl23|)). 

The operators of the quantum Brownian motion are introduced in appendix C, 
and satisfy the weak relations: 

dB\ dB t = ndt, dB t dB\ = (n + 1) dt, (13) 

dB t dB t = ndt, dB[ dB\ = (n + 1) dt, (14) 

and their tilde conjugates, with n being the Planck distribution function given in 
appendix B. (| and |) are the vacuum states representing the quantum Brownian 
motion. They are tilde invariant: (|~ = (|, |)~ = |). It is assumed that, at t = 
0, a relevant system starts to contact with the irrelevant system representing the 
stochastic process included in the martingale dM t . c 



3.3 Quantum Langevin Equation 

The dynamical quantity A(t) of the relevant system is defined by the operator 
in the Heisenberg representation: A(t) = VjT 1 ^) A Vf(t) where V^ _1 (i) satisfies 

dVf x {t) = Vf 1 ^) iHj t dt with Hj t dt = Hf, t dt + idM t dM t . 

In NETFD, the Heisenberg equation for A(t) within the Ito calculus is the 
quantum Langevin equation of the form 

dA(t) = i[H f (t)dt, A{t)] - d'M(t) [d'M(t), A(t)}, (15) 
with H f (t)dt = Vf 1 ^) Hf, t dt V f (t), and 

d'M(t) = VfHt) dMt V f (t). (16) 

Since A(t) is an arbitrary observable operator in the relevant system, ( [l5| ) can be 
the Ito's formula generalized to quantum systems. 

Applying the bra- vacuum ((1| = (|(1| to ( |l5| ) from the left, we obtain the 
Langevin equation for the bra- vector ((l|A(i) in the form 

d((l\A(t) = i{(l\[H s (t), A(t)]dt + ((l\A(t)n(t)dt - i((l\A(t) d'M(t). (17) 

In the derivation, use had been made of the properties (l|^4^(t) = (l|yl(i), (\d'W(t) = 
{\d'B(t), and ((l\d'M(t) = 0. 



3.4 Quantum Master Equation 

Taking the random average by applying the bra- vacuum (| of the irrelevant sub- 
system to the stochastic Liouville equation (|lTJ) , we can obtain the quantum master 
equation as 

(d/dt)\0(t))=-iH\0(t)), (18) 
with Hdt = (\Hf, t dt\) and |0(i)) = <|0 /(*)>■ 

c Within the formalism, the random force operators dBt and dB^ are assumed to commute 
with any relevant system operator A in the Schrodinger representation: [A, dBt] = [A, dB' t ] = 
for t > 0. 
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3.5 Stratonovich-Type Stochastic Equations 

By making use of the relation between the Ito and Stratonovich stochastic calcu- 
luses, we can rewrite the Ito stochastic Liouville equation ( |Tl| ) and the Ito Langevin 
equation ( p"5| ) into the Stratonovich ones, respectively, i.e., 

d\0 f (t))=-iH ftt dto\0 f (t)), (19) 

with H fit dt = H s dt + i(fldt + \dM t dM t ) + dM t , and 

dA(t) = i[H f {t)dt ° A(t)}, (20) 

with H f (t)dt = H s (t)dt + i{II(t)dt + \d'M{t)d'M{t)) + : d'M(t) : . 



3.6 Fluctuation- Dissipation Relation 

The fluctuation-dissipation theorem of the second kind for the multiple of martin- 
gales, dM t dMt, is determined by the criterion that there is no diffusive term comes 
out in the terms fldt + \dM t dM t appeared in Hf^dt in subsection |3.5| : 

dM t dM t = -2Tl D dt. (21) 

The origin of this criterion is attributed to the way how the Langevin equation 
was introduced in physics, as explained before, i.e., relaxation term and random 
force term were introduced in mechanical equation within the Stratonovich calculus. 
Therefore, there is no dissipative terms in stochastic equations of the Stratonovich 
type. We adopted this criterion in quantum cases. 



The operator relation (21) may be called a generalized fluctuation-dissipation 



theorem of the second kind, which should be interpreted within the weak relation. 



4 A System in the Rotating Wave Approximation 

4.1 Model 

We will apply the above formalism to the model of a harmonic oscillator embedded in 
an environment with temperature T. The Hamiltonian Hg of the relevant system is 
given by 77s = uia'a where a, and their tilde conjugates are stochastic operators 
of the relevant system satisfying the canonical commutation relation [a, a^] = 1, 
and [a, <r] = 1. The tilde and non-tilde operators are related with each other by the 
relation (l\a' = (l\d where (1| is the thermal bra-vacuum of the relevant system. 

Since we are interested in the system in the rotating wave approximation, we 
will confine ourselves to the case where the stochastic hat-Hamiltonian Ht is bi- 
linear in a, a^, dBt, dBj and their tilde conjugates, and is invariant under the 
phase transformation a — > ae l ° , and dBt — > dBt e 10 . This gives us the system of 
linear-dissipative coupling. 



Then, II r and II d consisting of 77 introduced in subsection 3.2 become 



IIr = + 7*7v). n D = 2n{n + v)^^ 1 , (22) 
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respectively, where we introduced a set of canonical stochastic operators 7^ = 
fia + va\ 7^ = at — a with fi + v = 1, which satisfy the commutation relation 
[7,,, 7+] = 1. The parameter v (or fi) is closely related to the ordering of operators 
when they are mapped to c-number function space with the help of the coherent 
state representation EU, i.e., v = 1 for the normal ordering, v = for the anti- 
normal ordering, and v — 1/2 for the Weyl ordering. The new operators 7^ and 7^ 
annihilate the relevant bra- vacuum: (1|7^ = 0, (1|7 = 0. 



4-2 Martingale Operator 

Let us adopt the martingale operator: 

: dM t : = : dM { ~ ] : + A : dM t (+) : (23) 

with : dM^ ] : = i(^dW t + j*dW t ) and : dAf t (+) : = -i{dW^ lv + dwf%). Here, 
the annihilation and the creation random force operators dW t and aW t are defined, 
respectively, by dW t = V2K(pdB t + vdB\), dWf = V2n(dB\ - dB t ). The latter 
annihilates the bra-vacuum (| of the irrelevant system: (|aW/ = 0, (\dW^ = 0. 
Note that the normal ordering : • • • : is defined with respect to 7's and dW's. 

The real parameter A measures the degree of non-commutativity among the 
martingale operators: [ : dM^ :, : dM^ : ] = —21Jjidt. In deriving this, we used 
the facts that 

dW t dW t = dW t dW t = 2k (n + v) dt, dW t dwf = dW t dW^ = 2ndt, (24) 

and that the other combinations are equal to zero. Note that [dWt, dWf] = 
2ndt sh,atdd be compared with @). There exist at least two physically attractive 
tfl, i.e., one is the case for A = giving non-Hermitian martingale: 



dM t = W2k 



(a f - a) dSl + t.c. , (25) 

and the other for A = 1 giving Hcrmitian martingale: 

dM t = iV2K[(aUB t - dB\a^ + t.c], (26) 

where t.c. stands for tilde conjugation. The former follows the characteristics of 
the classical stochastic Liouville equation where the stochastic distribution function 
satisfies the conservation of probability within the phase-space of a relevant system 
(see section ||) . Whereas the latter employed the characteristics of the Schrodinger 
equation where the norm of the stochastic wave function preserves itself—- In this 
case, the consistency with the structure of classical system is destroyed ErEJ. 
The fluctuation-dissipation theorem of the system is given by 

: dM t : : dM t := -2(Xn R + fl D )dt, (27) 

where we used the relations : dM 4 (_) : : dM t ( ~ } := -21J D dt, : dM 4 ( ~ } : : dM t (+) := 
-21J R dt and : dM ( t +) : : dM ( t +) :=: dM t (+) : : dM t ( ~ } := 0, which can be derived by 
making use of (EJ). 
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The hat-Hamiltonians of the model are given by 

H f<t dt= H s dt + i(l- X)fl R dt + dM t , (28) 

H] t dt = H s dt + i((l - 2\)h R - fl D )dt + dM t , (29) 

H f (t)dt = H s {t)dt + i(l - X)n R (t)dt+ : d'M{t) : . (30) 

4-3 Heisenberg Operators of the Quantum Brownian Motion 

The Heisenberg operators of the Quantum Brownian motion arc defined by 

B(t) = Vf\t) B t V f (t), B\t) = Vf\t) B\ V f (t), (31) 

and their tilde conjugates. Their derivatives dB*{t) = d{Vf l {t) Bf Vf{t)), (# 

: nul, dagger and/or tilde) with respect to time in the Ito calculus are given, 
respectively, by 

dB(t) = dB t + V2^ [(1 - A) v (a f (t) - a(t)) - Xa{t)] dt, (32) 

dB f (t) = dB\ - V2^ [(1 -X)n (a + (t) - a(t)) + Xa\t)] dt, (33) 
and their tilde conjugates. Then, we have 

dW(t) = dW t - X2K~/v(t)dt, dW l {t) = dW^ - 2n 1 l {t)dt. (34) 
Since, by making use of (|34|), we see that 

dM(t) = d'M(t) = i[^{t)dW t + ^{t)dW t } - iX[dW^ v {t) + dWf >(*)], (35) 

we know that the martingale operator in the Heisenberg representation keeps the 
property: (\dM(t)\) = 0. 

4-4 Quantum Langevin Equations 
The quantum Langevin equation is given by 
dA(t) = i[H s (t), A(t)]dt 

+K{(l-2A)( 7 *(*)[7*(t), A(t)]+f(t)[Ut), A(t)}) 

+[ 7 *(t), A(t)] lu (t) + lf(t), A(t)]%(t)}dt 
+2 K (n + is)[^(t), [ 7 *(t), A{t)]]dt 
-{[ 7 *(t), A(t)]dW t + [f(t), A(t)}dW t } 

+A{dW t *[ 7 „(t), A(t)]+dW?tf v (t), A(t)}} (36) 
= i[H s (t), A(t)]dt 

+^Ht)hAt), A(t)]+jHt)lUt), A(t)] 

+(1 - 2A)([ 7 *(t), A{t)] lv (t) + [f{t), A{t)Y lv (t)))dt 
+2 K (n + v)[y*{t), [ 7 *(t), A{t)]]dt 
-{[ 7 *(*), A(t)]dW(f) + [ 7 *(f), A(i)]cW(t)} 

+A{dW*(t)[ 7 „(t), A(t)]+rfly*(t)[ 7 ,,(t), (37) 
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) in ( p6|), and by means of that in the 
3) in l37|). 



with H s (t) = V f 1 {t)H s Vf(t) = H s (t) - Hs(t). Note that the Langevin equation 
is written by means of the quantum Brownian motion in the Schrodinger (the in- 
teraction) representation (the input field cl 
Heisenberg representation (the output field 

The Langevin equation for the bra- vector state, ((l\A(t), reduces to 

d((l\A(t) =i((l\[Hs(t), A(t)]dt 

- K {((l\[A(jk), at(t)]o(t) + ((l|ot(t)[o(t), A(t)]}dt 
+2/cn«l|[a(f), [A(t), a^(t)]]dt 

+ ((l\[A(t), a\t)]V2^ dB t + ((1\V2^ dBl[a(t), A(t)} (38) 
= i((l\[H s (t), A(t)]dt 

-K(l-2A){«l|L4(t), aHt)]a(t) + ((l\aHt)[a(t), A(t)]} dt 
+2«n«l|[a(t), [A(t), a^(t)]]dt 

+ ((l\[A(t), a\t)]V2ildB{t) + ({l\V2ildB\t)[a{t), A(t)}. (39) 

The relation between the expression ( J38| ) and ( ]39j ) can be interpreted as follows. 
Substituting the solution of the Heisenberg random force operators ( |32"| ) and ( |3^ ) 
for dB(t) and dB^(t), respectively, into (j39[), we obtain the quantum Langevin 
equation (Q) which does not depend on the non-commutativity parameter A. 
The Langevin equations for a(t) and a^(t) of the system reduce to 

da(t) = {-iuj - k) a(t)dt + dW t - 2(1 - \)vn [a\t) - a(t)] dt - \vdWf, (40) 
da\t) = {iu - k) a\t)dt + dW t + 2(1 - X)fXK [a\t) - 5(t)] dt + XfidW? . (41) 

Note that the last two terms in the above equations disappear when one applies 



to them. For A = 0, (|40|) and (J41|) become, respectively, to 

da(t) = -iuja(t)dt - nd\t)dt + dW t , 
— iuia) (t)dt — Ka{i)dt + dWt, 



where we put ji 



= 1/2, for simplicity. For A = 1, we get 

da(t) = —iuja(t)dt — K,a(t)dt + V2ndB t , 
do){t) (t)dt - na\t)dt + V^ndBj, 



(42) 
(43) 



(44) 
(45) 



which may correspond to (|2|). Applying = (1|(| to ( f40| ) and (J4l]) , we obtain, for 
any values of A, fi and v, the Langevin equations of the vectors {(1 |a(t) and ((l|a t (i) 
in the forms 



d((l\a(t) = -ioj((l\a(t)dt - n((l\a(t)dt - 
d((l\a^(t) = iuj{{l\a\t)dt - K{{l\a\t)dt 

Note that these have the same structure as those in ([ 
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2re«l|<m t , 
2n{{l\dB\. 



(46) 
(47) 



5 A System with Commutative Random Force Operators 

Let us investigate Boson system having x-X type interaction between environment, 
i.e., a system without the rotating wave approximation. The Hamiltonian of a 
harmonic oscillator can be written in the form H$ =■ p 2 /(2m) + moj 2 x 2 /2 : with 
x = wl/2mw(a + cv), p = — i^Jmuj /2{a — cv) where x and p satisfy the canonical 
commutation relation [x, p] = i. The normal ordering here, is taken with 

respect to a and cr. The relaxational and the diffusive parts in 77 = II R + II d are 
given, respectively, as 

il R = -in(x- x){p + p), II D = -2Kmui(n + 1/2) (x- x) 2 . (48) 

The martingale operator corresponding to x-X type interaction may have the 
form 

dM t = 2^Kmuj(xdX t - xdX t ), (49) 

with dX t = (dB t + <iB])/v / 2 where dBt and dB\ are the quantum Brownian motion 
defined in appendix C. Then, we have dX t dX t = dX t dX t = (n + 1/2) df which 
gives us the fluctuation-dissipation theorem 

dM t dM t = -2fl D dt. (50) 

The form of the martingale (^) was adopted by following the structure of micro- 
scopic interaction Hamiltonian of the x-X type. 

The stochastic hat-Hamiltonian Tif^dt for the stochastic Liouville equation (|1 
of the Ito type is given by 

Hf,tdt = H s + i{n R + n D )+dM t . (51) 

Then, the stochastic hat-Hamiltonian of the Stratonovich type becomes 

H f<t dt = H s dt + ih R dt + dM t , (52) 



where one does not see IId thanks to the fluctuation-dissipation theorem Q50| ) . We 
can also check that 

nj t dt = H s dt + i(II R - II D )dt + dM t . (53) 
The Langevin equation has the forms 

dxit) = —p(t)dt + n(x(t) - x(t))dt, (54) 
m 

dp(t) = -muo 2 x(t)dt - n(p(t) + p(t))dt 

+4:iKmu}(n + l/2)(x(t) - x(t))dt - 2^ nmw dX t . (55) 

Applying (1| to ( ^) and (|55|), we have the Langevin equation for (l\x(t) and 
(l|p(i) in the forms 

d(l\x(t) = m- 1 (l\p(t)dt, (56) 
d(l\p{t) = -muj 2 (l\x{t)dt - 2n{l\p(t)dt - 2y/^U (l\dX t , (57) 

10 



respectively. They can be written in terms of a and a' as 

d(l\a(t) = -iu(l\a(t)dt - k{1\ [a(t) - a f (f)] dt - (dB t + dB^ , (58) 

d(l| a t( t ) = iu(\\a){t)dt - k(1\ [a f (t) - a(t)] dt - iy^(l\ (dBj + dB t ^j . (59) 

If we take the rotating wave approximation at this stage the coefficients in front of 
the quantum Brownian motion are not equal to those appeared in ( [44| ) and (p5[). 
It may indicate that a naive procedure of taking the rotating wave approximation 
will not give us correct results. It might also be related to the renormalization 
procedure needed to derive stochastic differential equations Jor the system with 
x-X type interaction from a microscopic Heisenberg equation Ej. 



6 Concluding Remarks 

We have revealed that the non-commutativity among dM^ and dM^ appeared 
in the martingale operator of the model within the rotating wave approximation 
affect the relaxation part of the stochastic hat-Hamiltonian. When the measure 
A of the non-commutativity has the value A = 1, the hat-Hamiltonian becomes 
Hermite, and therefore, it looks like being related to a microscopic description. 
On the other hand, for A = 0, the system of the quantum stochastic differential 
equations has the same structure as that of classical mechanics, and it is related to 
a semi-macroscopic description. As has been shown in this paper, Hermiticy of the 
hat-Hamiltonian is realized thanks to the non-commutativity between dM^ ^ and 
dMi +) . 

Does this mean that the system with commutative martingale does not have any 
microscopic realization? As an example of system with commutative martingale, we 
studied the system corresponding to the quantum Kramers equation which has x-X 
type interaction Hamiltonian between the relevant system and environment system. 
Since, there is no non-commutative parts in the martingale operator, this system 
cannot have an Hermitian stochastic hat-Hamiltonian. The non-commutative parts 
appears when one takes the rotating wave approximation to the interaction Hamil- 
tonian. Does dissipation originate in the approximation causing non-commutative 
character in martingale? Can quantumness, appeared in this way, be the origin of 
dissipation? On the contrary, the following question arises naturally. Is it always 
possible to put all random force operators to be commutative? 

There are still a lot of problems to be resolved before we know the origin of 
dissipation. However, with the help of NETFD, we can see the problems from a 
unified viewpoint which may provide us with good prospects for further develop- 
ments. Introducing the parameter A in the martingale term as given by (|23|), we 
can transform the equation to the non-Hermitian version by shifting A — > (see 
). In other words, it seems that the non-commutativity is renormalized into the 
relaxational and diffusive terms. 

Substituting the solution of the random force operators (^) and ( |33|) in the 
Heisenberg representation (the output field) into (|37j), we have the Langevin equa- 
tion ( |3rj| ) expressed by means of those in the Schrodinger (or, more properly, the 
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interaction) representation (the input field). Note that the Langevin equation ( j3q ) 
for the bra-vector state ((l\A(t) does not depend on A when it is represented by the 
random force operator in the Schrodinger representation (the input field). 

We are intensively investigating what is the physical meaning of the renormal- 
ization of non-commutativity by changing the parameter A. The relation between 
the present argument and the procedure of the coarse graining is under investiga- 
tion. Related to the system with commutative random force operators, a micro- 
scopic derivation of quantum stochastic equations corresponding to the quantum 
Kramers equation are in progress E3. There, an appropriate renormalization is re- 
quired in accordance with the separation of two time-scales, i.e., microscopic and 
macroscopic time-scales. Without the renormalization, one gets quantum stochas- 
tic equations in the rotating wave approximation, which do not correspond to the 
system described by the Kramers equation. Including these studies, the further 
progress will be reported elsewhere. 
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Appendix A 

The definitions of the Ito and the Stratonovich stochastic products are given, 
respectively, by 

X t alY t = X t (Y t+dt -Y t ), dX t Y t = (X t+dt - X t )Y t (60) 

and 

X t odY t = Xt+dt + Xt {Y t+dt -Y t ), dX t oY t = (X t+dt -X t ) Yt+dt 2 +Yt (61) 

for arbitrary stochastic operators X t and Y t . From ([30|) and (^lj), we have the 
formulae which connect the Ito and the Stratonovich products in the differential 
form 

X t o dY t = X t dY t + (l/2)dX t dY t , dX t o Y t = dX t Y t + (l/2)dX t dY t . (62) 
Appendix B 

The time-evolution of the thermal vacuum |0(t)), satisfying the quantum master 
equation (081) with the hat-Hamiltonian for the semi-free system specified by H$ = 
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uja^a and (|22|), is given by 

\0(t)) =cxp{[n(<)-n(0)] 7 *7*}|0) ; (63) 

where the one-particle distribution function, n{t) — ((l|at(i)a(i)|0)), satisfies the 
kinetic (Boltzmann) equation of the model: dn(t)/dt = — 2k [n(t) — n) with the 

Planck distribution function n — (e u / T — l) 1 . Here, T is the temperature of 
environment system. 

Appendix C 

Let us introduce the annihilation and creation operators bt, b\ and their tilde con- 
jugates satisfying the canonical commutation relation: 

[bt, b\,] = 5(t - t'), [b t , %] = S(t - t 1 ). (64) 

The vacuums (0| and |0) are defined by b t \0) = 0, 6 t |0) = and (0\b\ = 0, (0|& t f = 
0. The subscript or the argument t represents time. 

Introducing the operators B t = J * dt dB t > — C dt' b v , B\ = J*' dt dB\, = 
L dt' b\, and their tilde conjugates for t > 0, we see that they satisfy B(0) = 0, 
B^(0) = 0, [B s , B\] = min(s, t), and their tilde conjugates, and that they annihilate 
the vacuums |0) and (0|: dB t \0) = 0, dB t \0) = 0, (0|dB t f = 0, {0\dB f t = 0. These 
operators represent the quantum Brownian motion. 

Let us introduce a set of new operators by the relation dC^ = B^dB" with 
the Bogoliubov transformation defined by 

8""= f 1 ^ - ") , (65) 
where n is the Planck distribution function. We introduced the thermal doublet: 
dB? =1 = dB u dB^ 2 = dB$, dBT 1 = d B l, dB~t =2 = ~ dB t, (66) 

and the similar doublet notations for dC t M and dC^. The new operators annihilate 
the new vacuum (| and |): dC t \) = 0, dC t \) = 0, (\dCj = 0, {\d(j\ = 0. 

We will use the representation space constructed on the vacuums (| and |). 
Then, we have, for example, 

(|dBt|) = {\dB\\) = 0, (\dB\dB t \) = ndt, (\dB t dB}\) = (ft + 1) dt. (67) 
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